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Aims. We study jet-driving protostars and T Tau stars to characterize their X-ray emission. We seek soft spectral components that may be due 
to shock emission, and search for soft, shock-induced emission displaced from the stellar position. 

Methods. We study two stellar samples, the first consisting of lightly absorbed T Tau stars with strong jets, the other containing protostars with 
disks seen nearly edge-on. The former sample was observed in the XMM-Newton Extended Survey of the Taurus Molecular Cloud (XEST), 
. while Chandra archival data provided observations of the latter. 
<^^j Results. We confirm the previously identified peculiar spectrum of DG Tau A and find similar phenomenology in GV Tau and DP Tau, sug- 
gesting a new class of X-ray spectra. These consist of a lightly absorbed, very soft component and a strongly absorbed, very hard component. 
■ The latter is flaring while little variability is detected in the former. The absorption of the harder component is about an order of magnitude 
~"^h higher than expected from the optical extinction assuming a standard gas-to-dust mass ratio. For the absorbed protostars, only the hard, stellar 
£s . X-ray component is found. 

I ' Conclusions. The flaring hard component represents active coronal emission. Its strong absorption is attributed to mass inflow from the ac- 
O cretion disk. The optical extinction is small because the dust has sublimated at larger distances. The weakly absorbed soft component cannot 
4— i originate from the same location. Because the stars drive strong jets, we propose that the X-rays are generated in shocks in the jets. We find 
^ that for the three peculiar X-ray sources, the luminosity of the soft component roughly scales with the equivalent width of the [O I] A 6300 line 
. . formed in the jets, and with the mass outflow rate. In the more strongly obscured protostars, the soft component is entirely absorbed, and only 
^ ■ the hard, coronal component penetrates the envelope or the near-edge-on disk. 
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1. Introduction 

Pre-main sequence stars show various signs of accretion and 
outflow, such as circumste llar disks detected at radio and mil- 
limeter wavelengths (e.g., ISimon et al.ll2000|). molecular o ut- 
flows observed in molecular lines (e.g., Bachiller 1996 for 
a review), and accomp anying optical (e.g.. iHirth et all 1997: 
lEisloffel & Mundtll 19981) and radio jets (e.g.. lAngladalll995l 
The most visible manifestation of outflows are the optically 
visible jets that may be excited close to the star, in internal 
shocks along the mass stream, or at distances up to several 
parsecs as the fast stream encounters the interstellar medium 
where it shock-ionizes the gas in Herbig-Har o objects (fo r 
a review of Herbig-Haro flows, see Reipurth & Bailvl l200ll) . 
Under ideal circumstances (low extinction, strong ionization), 



jets can be identifi e d at distances as close as 0.1" to the 
star JBacciotti et all 120021) . The same compact jets are also 
routinely detected at radio wavelengths, where the emission 
mechanism is thought to be b remsstrahlung from the shock- 
heated gas jRodrfguezl 1 19951 lAngladaL Il995h . Radio bright- 
ness temperatures suggest overall gas temperatures of order 
10 4 K. This picture is ambiguous, however, as a number of 



tion or synchrotron-like spectral shapes 


Yusef-Zadeh et al. 


199fi ICuriel et al. 


119931: iHughesl 11997 


iRav et all Il997 


Smith et all 120031 


Loinardet al., 2005]). Magnetic fields may 
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thus play a role not only in the launching of jets, but in their 
propagation as well. 

Two principal families of models for the jet-launching 
mechanism ha ve been proposed: The magnetocentrifu- 
gal mechanism ( B landford & Pavnei Il982t iKonigl & Pudritzl 
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2000) posits that outflows are launched on the disk along 
outward-bent, open magnetic fields owing to the increasing 
centrifugal forces along the field lines that are co-rotating with 
their d isk footpoints. In contrast, the X-wind model f Sh u et all 
l2000l) assumes that outflows are launched from magnetized re- 
gions between the star and the disk. Magnetohydrodynamic 
simulations show that magnetic reconnection and subse- 
quen t plasmoid ejec tion may lead to field-aligned outflows 
llHavashiet al.Lll996h . 

Accretion and outflow processes are prone to producing 
X-rays, given that shocks with shock jump velocities of or- 
der several hundred km s _! are possible and sometimes ob- 
served. Accretion onto the star essentially extracts gravitational 
energy and transforms it into thermal energy through shocks . 
The relevant theory has bee n presented by lUlrichl dl976l) . 
ICalvet & Gullbringl(ll998l) and lLamzinN 19991) . Material falling 
along magnetic fields toward the stellar surface with nearly 
free-fall velocity, Vff, develops strong shocks near or within 
the photosphere, heating plasma to temperatures of T m 
3v\nm p l\6k x 5.2 x lO^M/R [K] where M and R are the 
stellar mass and the stellar radius, respectively, in solar units, 
/1 is the mean molecular weight, m p is the proton mass, and 
k is the Boltzmann constant. For T Tau stars, heating to a 
few MK is thus en tirely plausible (as, e.g., proposed for TW 
Hva. iKastner et al1l2002l) . but the temperatures in excess of 
10 M K that often dominate X-ray emitting pla smas in T Tau 
stars JSkinner et alll2003llPreibisch et alll2005l) cannot be due 
to accretion shocks. Rather, high electron temperatures in pre- 
main sequence stars are conventionally attributed to magne- 
tized plasma trapped in corona-like or large-scale "magneto- 
spheric" magnetic fields around the stars. 

Shocks also form in jets and outflows, in particular in 
Herbig-Haro (HH) objects. T he relevant theory and a simple 
model have been discussed by Raga et all J2002I) . The strong- 
shock temperature can be expressed as T =s 1.5 x 10 5 Vj ()() K 
(for fully ionized gas) where vioo is the shock speed relative to 
a target in units of 100 k m s" 1 . Jet speeds are typically of or- 
der v = 300 -500 km s " 1 feisloffel & Mundtlll998tfAngladal 
ll995UBallve7aill2003l) . allowing for shock speeds of similar 
magnitude. If a flow shocks a standing medium at 400 km s , 
then T « 2. 4 MK. 

Observations have explicitly demonstrated that faint, 



soft 


X-ray e 


mission originates from some HH objects 


(Pravdo et al.,2001, 


2004; Favata et al., 2002; Tsuiimoto et al., 


2004: 


; Pravdo & Tsuboi. 2005i|Grosso e{ al] 


2006) 


.Bally et alJ 



(2003) used a Chandra observation to show that X-rays form 
within an arcsecond of the protostar L155 1 IRS-5 while the star 
itself is too heavily obscured to be detected. As this example il- 
lustrates, the jet-launching region of powerful protostellar jets 
is usually inaccessible to optical, near-infrared or X-ray stud- 
ies of protostellar outflow sources due to excessive absorption. 
However, there are a number of optically-revealed , classical T 
Tau s tars that drive appreciable jets and outflows dHirth et all 
Il997l) . Their jets were initially dubbed "micro-jets" although 
recent studies have shown them to reach parsec scales and the 
mass loss rates com pete with those of more embedded sources 
jMcGroartv & RavihofrA . 



One of the most prominent examples, the classical T Tau 
star DG Tau A, has recently obtained scrut iny with high- 
resolution Chandra X-ray observations. IGiidel et all {2005) 
discovered a new type of X-ray spectrum in which a lightly ab- 
sorbed spectral component from a rather cool (2-3 MK) plasma 
is complemented by a strongly absorbed component from a 
very hot (20-30 MK) plasma. The latter was interpreted as 
being related to the usual coronal/magnetospheric source of- 
ten seen in T Tau stars, subject to some excess absorption. In 
the light of this interpretation, the cooler plasma is unlikely to 
be located close to the magnetospheric component because the 
source of absorption would likely affect it as well. The strong 
jets and high mass-loss rates of DG Tau A motivated the sug- 
gestion that the soft X-rays are produced in shocks in the ac- 
celeration region of the jet (the "base" of the jet), which is un- 
resolved in Chandra or XMM-Newton images. Complementary 
evidence was seen in faint and very soft X-ray emission along 
the optical jet, outside the stellar point-spread function and co- 
spatial with a strong bow-shock feature seen in the optical. The 
spectral distribution of the counts in the jet source was reminis- 
cent of the count distribution in the soft stellar component. A 
spectrum requiring two absorbers was also reported for a pro- 
plyd in the the Orion Nebula Cluster, wit h a variable ha r d com - 
ponent but a consta nt soft component ( Kastne r et all 2005). 
ISkinner et all (120061) have detected a similar two-component 
spectrum in the strongly accreting FU Ori, the prototype of 
FU Ori-type stars. The cooler, weakly absorbed component, 
however, reached temperatures more common to T Tau stars 
(» 8 MK) and was therefore identified with a coronal com- 
ponent, while the very hot component was again strongly ab- 
sorbed. 

The recently conducted XMM-Newton Extend ed Survey of 
the T aurus Molecular Cloud (XEST henceforth; IGiidel etall 
2006) offers unique access to numerous stellar objects in the 
Taurus Molecular Cloud (TMC). It covers most stellar concen- 
trations in the TMC molecular filaments and offers high-quality 
spectra and light curves for any "typical" T Tau star in the field 
of view. XEST also covers a number of well-studied jet-driving 
T Tau stars and protostars. The observations confirm the ini- 
tial discovery of the anomalous DG Tau A and add substantial 
evidence for a jet-related origin. We find further clear exam- 
ples of these two-absorber X-ray sources (TAX sources hence- 
forth). The present paper describes these new observations and 
presents a detailed discussion on the possible origin of the X- 
ray emission. For a comparison, we also present X-ray spectra 
of jet-driving sources that are deeply absorbed by their thick 
circumstellar disks seen nearly edge-on. The soft component, 
if present, would thus be absorbed unless it is displaced suffi- 
ciently from the star to make it separate from the stellar X-ray 
image. These complementary data are based on high-resolution 
Chandra observations. 

The plan of the paper is as follows. In Sect.|2] we briefly de- 
scribe our data analysis and introduce the targets of this study. 
Sec.|5]presents the results that are further discussed in Sect.|4] 
We conclude in Sect. |5] 
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Parameter 



XEST number 
ObsID (XMM-Newton) 
Start time (UT) 
End time (UT) 
Exposure time (s) 
Boresight RA (J2000.0) 
Boresight 6 (J2000.0) 
Star: RA (J2000.0)" 
Star: 6 (J2000.0) 



DG Tau A 



02-022 
0203540201 
2004-08-17 06:08:10 
2004-08-17 17:32:46 
41076 

04h 27m 19.6s 
26deg 09' 25" 
04h 27m 04.70s 
26deg 06' 16'.' 3 



GV Tau 



13-004 
0203541301 
2004-08-25 11:11:23 
2004-08-25 19:59:19 
31676 

04h 29m 52.0s 
24deg 36' 47" 
04h29m 23.73s 
24deg 33' 0(X'3 



DP Tau 



10-045 
0203542201 
2005-03-05 05:56:38 
2005-03-05 14:44:30 
31672 

04h 42m 20.9s 
25 deg 20' 35" 
04h 42m 37.70s 
25 deg 15' 37'.'5 



CWTau 



20-046 

0203542001 

2004-09-12 07:04:43 

2004-09-12 15:52:37 

31674 

04 14 12.9 

28 12 12 

04h 14m 17.00s 

28 deg 10' 57'.'8 



Parameter DG Tau B FS Tau B HH 34 IRS HH 111 IRS 

XEST number C3-2 11-054 = C2-1 

ObsID ( Chandra) 4487 4488 4489 4490 

Start time (UT) 2004-01-11 02:58:51 2003-11-08 12:57:58 2003-12-05 08:59:23 2004-11-21 17:05:50 

Endtime(UT) 2004-01-11 11:52:21 2003-11-08 21:56:34 2003-12-05 17:38:31 2004-11-22 01:36:42 

Exposure time (s) 29717 29674 28590 28795 

Boresight RA (J2000.0) 04h 27m 02.3s 04h 22m 00.2s 05h 35m 32.07s 05h 51m 54.97s 

Boresight S (J2000.0) 26 deg 04' 56" 26 deg 58' 07" -06 deg 26' 50r0 02 deg 48' 55'.'3 

Star: RA (J2000.0)* 04h 27m 02.55s 04h 22m 00.70s 05h 35m 29.84s 05h 51m 46.31s 

Star: S (J2000.0) 6 26 deg 05' 3079 26 deg 57' 32'.' 50 -06 deg 26' 58'.' 40 02 deg 48' 29T 72 

a For origin of coordinates, see Giidel et alH2006l) 
b For origin of coordinates, see Tables|5|and|6| 



2. Targets and observations 

2.1. Data reduction 

The observations reported here were obtained as part of the 
XEST project and have been subj ect to a data reduction proce- 
dure that is descri bed in detail inlGtidel et alJ (120061) . In brief, 
the XMM-Newton Jjansen et all 1200 ll) data have been obtained 
using all three EPIC cameras (MOS and PN: lTurner et alJ200ll 
IStriider et alJl200ll) . with the medium filter inserted. We ap- 
plied extraction regions that maximize the signal-to-noise ra- 
tio; their radii amount to 22", 15", 13", and 12" for our tar- 
gets DG Tau A, GV Tau, DP Tau, and CW Tau, respectively. 
Standard response matrices distributed by the XMM-Newton 
project were used, and ancillary response files were created for 
each source. An observing log is given in Tabled For DG Tau 
A and GV Tau, we also extracted the simultaneous light curves 
from the Optical Monitor (OM; Mas on et alJl20 01). The OM 
data reduction is described in Audard et al. ( 2006). 

For spectral ana lysis, we used the XSPEC software pack- 
age dArnaudll 19961) . We applied the vapec thermal collisional 
ionization equilibrium models that include emission lines and 
bremsstrahlung continua, combined with the photoelectric ab- 
sorption model, wabs, based on absorption cross sections by 
iMorrison & McCammonldl983l) . 

The Chandra exposures described here were reduced us- 
ing standard strategies for ACIS. The observations of DG Tau 
AB, FS Tau AB (= Haro 6-5 AB), HH 34, and HH 1 1 1 used 
the "Very Faint Mode" of ACIS to minimize background con- 
tributions, while the GV Tau (= Haro 6-10 AB) observation 
used the "Faint Mode". The data were reduced in CIAO vers. 



3.0.2 following the standard analysis threads 1 . These proce- 
dures included corrections for charge transfer inefficiency and 
afterglow, and selection of good time intervals. Spectra were 
extracted from level 2 event files with specextract which also 
produces response matrices and ancillary response files. 

2.2. Jet-driving T Tau stars 

The XEST survey contains several jet-driving , optically re- 
vealed classical T Tau stars jHirth et all Il997l) . in particular 
DG Tau A, GV Tau AB = Haro 6-10 AB, DP Tau, CW Tau, 
XZ Tau, UZ Tau, and DD Tau. One further object, HN Tau, 
was observed with th e Chandra H RC-I but was recorded only 
as a very faint source (iGtidel et alll2006l) . 

At least three XEST targets, namely DG Tau A, GV Tau 
A, DP Tau, and possibly a fourth, CW Tau, reveal the same 
anomalous spec tral phenom e nolo gy previously described for 
the first object JGtidel et all l2005h . Further, UZ Tau reveals 
a peculiarly soft spectrum, while DD Tau shows a very hot 
source that is subject to an anomalously high photoelectric ab- 
sorption, jud ged from the visual extinction which is relatively 
small (Giidel et al., 2006). The present paper concentrates on a 
discussion of the first three unambiguous cases. 

Basic properties of the four peculiar X-ray sources are col- 
lected in Table 13 The table lists spectral type, infrared clas- 
sification ("YSO class"), T Tauri type according to the Ha 
equivalent width, t he stellar effective temperature T e e, the stel- 
lar mass based on Sies s et alJ (|2000) evolutionary tracks, the 

1 http://cxc.harvard.edu/ciao/guides/ 
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Table 2. Fundamental parameters of optically revealed stars 



Parameter 


DG Tau A 


Ref" 


GV Tau A 


Ref 


GV Tau B 


Ref 


DP Tau 


Ref 


CWTau 


Ref 


Spectral type 


K5-M0 


4, 5, 7, 8, 15 


K3-7 


4, 8 






M0-0.5 


5, 8 


K3 


7 


YSO class 


II 


4, 8 


I 


4, 6, 8 


I 




II 


4, 6, 8, 18 


II 


4, 6, 8 




FS* 


6, 18 


















T Tau type 


CTTS 


5, 17 


CTTS 


8 






CTTS 


5 


CTTS 


17 


r eff (K) 


4205 


7 


4730 


7 


5150 


12 


3778 


14 


4730 


7 




4775 


8 


4000 


8, 12 






3800 


8 






Mass (M Q ) 


0.9-1.8 


7, 8 r 


0.7-2.1 


4, 8 C 


2 


12 


0.52-0.53 


4, 8'' 


1.4 


4 C 




0.67 


15 




















0.88 


5 


















Radius (R B ) 


2.5 d 




2.8 rf 


... 


3.1 


12 


l.05 d 




L6 d 




W([0 I]) (A) 


11-22 


15 


3.8 


2 






0.7 


2 


3.5 


2 




















7.9 


15 


s -'* V-MD/ 


3.62 


8 


1.8 


8 






0.2 


14 


1.1 


7 




1.7 


7, 15 


1.3 


12 






0.3 


5 


0.76 


5 




1.15 


5 


















/ ,U„1 I / J/~\ 1 

A -'bOl V*- 1 ©/ 


6.4 

8 


4 
1 


6.98 


4 


6.1 


12 


0.7 


4, 8 


2.7 


4 


Av (mas) 


1.41 


5 


3.3 


11 


49 


12 


1.26 


5 


3.4 


15 




2.20 


10 


5.6 


12 






1.46 


4 


2.29 


4 




3.32 


8 


12.1 


8 


35 


13 


6.31 


8 








3.2 


15 


















^xj filing _/ 


0.36 


7 










0.41 


14 


0.55 


7 


losM,,.,. fMr»vr -1 ) 


-7.34 


5 


-6.71 


8 






-8.50 


5 


-7.99 


5 




-6.13 


8 










-6.92 


8 








-5.7 


15 


















lo2M,.„t (Mnvr -1 ") 

out \ irj -QJ L / 


-6.19 


8 


-6.54 


8 






-7.42 


8 


-7.1 


15 




-6.5 


15 




















-6.1 


10 




















-6.6 


16 


















Mawv (Mm) 


0.02 


6 


0.003 


6 






<0.0005 


6 


0.002 


6 


age (Myr) 


1.3-2.2 


7, 8 C 


0.89-0.98 


4, 8 C 






1.5-7.2 


8, 14 c 


7 


4 C 


Vrad (km S -1 ) 


250 


2 


80-300 


2 






90-110 


2 


110 


2 


v pm (kms _1 ) 


194 


3 


















A Vrad (km s -1 ) 






90 


2 






86 


2 






Pm (d) 


6.3 


19 


< 5.6 e 








< 2.8'' 




< 2.9'' 




v sin i (km s ) 


28.6 


8 


25.3 


8 






19.2 


8 


27.4 


8 



References: llCohen & Kuhill 19791) : 2 lHirfh et alJ J 19971) , radial velocities (along line of sight) for the hig h-velocit y outflow component; 
3 Doue ados et alJ 1 2000). proper motion for DG Tau A: 4lKenvon & HartmanrHl995[): 5lWhite & GhezH200lh: 6lAndrews & W illiams 
| 2005|) : 7 lBriceno et alJ <2002l): 8lwhite & Hillenbrand l2004 ; 9 2MASS.ICutri et alJ <2003l): lOlMuzerolle et alJ <1998l): 1 llLeinert et al 
120011): imoresko et alJ ll997h: 13lMenard et alJ <1993h: 14lLuhmarj 12004 : 15 lHartigan et alJ Jl995h : 16 lBacciotti et alJ <2002h : 17 
lKenvonetalJ<1998l) : 18 lHartmann et alJl2005h : 19 lBouvier et alj<1993l) 
FS = flat-spectrum source (also I/II) 



Ages and masses calculated based on parameters given in the cited references, usinglsi- 
Radius calculated from T sB and L bol 

Upper limit to rotation period calculated from radius and v sin i 



; et alJfcOOd) evolutionary tracks 



stellar radius computed from the stellar luminosity and T^, 
the equivalent width of [O il /i6300 lines (that can be formed 
in jets, Hirth et al. 1997), the stellar photospheric luminos- 
ity L„, the bolometric luminosity Lbol (integrated from the in- 
frared+optical spectral energy distribution, including contribu- 
tions from disks and envelopes), the extinctions Ay and Aj, 
the mass accretion rate M acc , the mass outflow rate M out , es- 
timates of the total disk mas s Mdisk, the estimated age as de- 
rived from lSiess et al 1 ll2000l) evolutionary tracks, the jet radial 
velocity v ra d and its dispersion Av, a d, the rotation period (de- 



rived from v sin i), and the projected equatorial velocity v sin 2. 
For DG Tau, we also list a value for the jet proper motion v pm . 
Multiple values are given if they are significantly different, as 
reported in the literature. A few notes on the individual stars 
follow. 

DG Tau A is a classical T Tau star whose infrared spec- 
trum, h owever, shows a rare, flat infr ared spectral energy distri- 
bution (Andrews & \yilli amsLl2005l) . It drives a very energetic 
jet (|Mundt & Fried, 1983) simila r to Class I prot ostars, includ- 
ing a counter-jet jLavallev-Fouauet et al.l l2000). The jets can 
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be traced out to a distance of 14.4' dMcGroartv & Ravi Eo04). 
DG Tau is one of the most active CTTS known, and ranks 
among the CTTS with the highest mass accretion and mass out- 
flow rates, competing with well-studie d protostellar jet sources 
dHartigan et al.Lll995t iBacciotti et all l2002l) . The star has in- 
deed been sugg ested to be a tran sition object between proto- 
stars and CTTS (P vo etall E()03). On the other hand, the stel- 



lar visual extinction is rather modest. 


with reported Ay in the 


range of 1.4-3.3 mag (Hartigan et al., 


1995; Muzerolle et al., 


1998; White & Ghez, 2001; 


White & Hillenbrand, 2004). A 



much lower visual extinction has been determined for the jet's 
Herbi g-Haro objects, namely Ay = 0.39 mag dCohen & Fulleil 
1985)? 

The jet is also known from radio observations 

dCohen & Biegingl Il986l) . High-resolution studies have 
shown that a narrow jet ploughs through progressively slower 
and wider outflow structures. A few arcseconds from the star, 
the proper motion of the jet indicates velocities in the plane o f 
the sky of 194 + 20 km s _] , on avera ge ([ Dougados etaim OOO). 
with a maximum of > 360 km s~ l (Eisloffel & Mundt, 1998). 
The jet radial velocity is 210-250 km s _I , compatible with 



the reported inclination angle of 38 deg ( Eisloffel & Mundt 



19981) . Bulk gas speeds rea ch 500 km s 1 ( B acciotti et al 
2000[ iBeristain et all 1200 lh. w ith FWHM line widths of 
100-200 km s~' dPvo et alll2003l) . 

The high-velocity jet shows bow-shock like structures 
out to distances of s everal arcsecs (also known as HH 158, 
iMundt & Friedlll983l) bu t can be followed down to 0.1" from 
the star dBacciotti et ail l2000t IPvo et all l2003h . Recently, 



tne star (tfacciotti et i 
IBacciotti et alL 2002 ) 



measured rotation of the jet around 
its flow axis, which has subsequently been used to infer 
the origin of the jet in the inner disk (0.3-4 AU for the 
lower-velocity component, and possibly in the X-wind re- 
gion _a^_die_jnner^isk^ge_a^_pT AU_for the high-velocity 
jet; lAnderson et aill2003t IPvo et alll2003l) . The jet mass-loss 
rate has been estimate d by various authors in the range of 
(2.4-7.9)xl0- 7 M o yr- 1 dHartigan et allll995tlMuzerolle et all 
ll998UBacciotti et alll2002HWhite & Hillenbrandll2004 . This 
rate amounts to ab out 30% of the accret i on rate, M acc = (0.8 - 
2) X l Q- 6 M^vr-' JHartigan etallll995l IWhite & HiUenb rand 
120041 IWhite & Ghe7»200ll give a somewhat lower accretion 
rate). Sh ocks seem to be the princ ipal source of heating and ex- 
citation dLavallev-Fouauet et all Eoool) . We note that DG Tau 
B, a jet-driving protostar, is located at a considerable distance 
from DG Tau A (separation « 50"), much larger than the PSF 
of the EPIC cameras. We will discuss this second source further 
below. 

GV Tau = Haro 6-10 AB is a close, embedded binary sys- 
tem with a separation of 1.3", consisting of an optically re- 
vealed T Tau-like star and a deeply embedded "infrared com- 
panion" considered to be as a protostar. The infrared spectrum 
of the system is that of a protostar (Class I, Tabled). This bi- 
nary has received part icularly detailed study in the radio range 
( Rei purth et all 120041). The jet shows r adial velocities in the 
range 80-200 km s" 1 dHirth et all f 19971) . 

DP Tau drives a jet that has been recorded to a distance 
of at least 27", and possibly to 111 ". The jet is oriented rela- 
tively close to the plane of the sky (Mundt & Eisloffel, 1998) 



although radial velocities of the jet still reach about 1 10 km s 1 
(Hirth et all 1 19971 ). The st ar shows strong [OI] emission from 
the iets dHirth et allfl997i) . 

CW Tau drives a bipolar outflow that can b e traced out to 
at least 4-6" from the star dHirth et all Il994h. with Her big- 
Haro objects found out to 16.8' (Mc Groartv & Ravi EooH) . A 
dedicate d study of this jet was presented by the latter authors. 
iDougados et all (EoOO) reported a highly collimated, transver- 
sally resolved jet structure. 

2.3. Jet-driving protostars 

We took advantage of the high spatial resolution of Chandra 
to obtain images of four protostellar sources that drive strong, 
collimated jets and outflows. Their circumstellar disks are seen 
nearly edge-on so that any soft emission will be strongly ab- 
sorbed unless it is originating sufficiently far away from the 
star. This allows us to distinguish between X-ray sources co- 
incident with the protostars and so fter X-rays from internal 
shocks in the jets (Ball v et all 120031) and will thus assist in our 
interpretation of the soft components seen in the T Tau sources. 
This geometry also permits us to probe the disks and envelopes 
themselves, by studying photoelectric absorption. 

Two objects, DG Tau B and FS Tau B = Haro 6-5B, are part 
of the XEST survey, but we also discuss here two further stars 
that are located in the Orion region at larger distances, namely 
the protostellar sources of the HH 34 and the HH 111 jets. 

We used the wavdetect algorithm included in the Chandra 
CIAO software to identify relevant X-ray sources and to de- 
termine their positions. The algorithm was applied to images 
containing photons in the 0.3-10 keV range, the inspected field 
containing approximately 10 6 standard ACIS pixels of size 
0.5" x 0.5". We set the sigthresh parameter at a value of 10~ 6 , 
which means that only one X-ray detection in the field is likely 
to be spurious. 

To optimize the positional information for the X-ray 
sources, we registered, a s far as pos s ible, f ield sources against 
2MASS catalog entries dCutri et all Eo03) but found that the 
positional agreement is extremely good, with typical offsets 2 
of < 0.4", which is not significant. We provide some specific 
information on the four sources below. 

DG Tau B is an embedded, envelope-dominated, 
jet-dr iving protostar d escribe d in detail b y, amon g oth- 
ers. IMundt & Fried] d 1983b . IMundt et all d 19871) . and 
Eisloffel & Mundd (1998). Near-infrared imaging shows 
two narrow cavities and a thic k dust-lane seen nea rly edge-on, 
but probably not the protostar ( Padg ett et all 19991 = P99). DG 
Tau B is considered to be in a Class I stag e jwalSQn_e_Lall 
120041). It was detected as a radio source by Rod riguez et all 

FS Tau B = Haro 6-5B is a jet-driving embedded source 
about 20" west of the classical T Ta u binary system FS Tau 
A. It was discovered as a jet source by Mundt et al. ( 1984) and 
describ ed i n detail by Mundt et a[| dl987h . lEisloffel & Mundl 
(1998), and lKrist et al. (1998). Near-infrared imaging shows it 



2 The 90% source location error circle in Chandra has a radius of 
about 0.5", see Chandra Proposes' Observatory Guide v.8. 
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as a disk-dominated protostar in which a thick dust lane, seen 
nearly edge-on, separates two nebulae considered to be the il- 
luminated surfaces of the flared disk, with a widely open cavity 
(P99). FS Tau B is probably in a rather evolved stage of Class I 
evolution, in transition to becoming a TTS. In the near-infrared, 
the star itself seems to be detecte d (P99). FS Tau B was found 
as a radio source bv lBrown et al.l |Tj?85). 

The HH 34 s ource drives one of the finest, highly colli- 
matedjets (Biihrke et allll988l) . The optical jet can be followed 
to within 1" of the Cl ass I protostar . The s tar itself has been de- 
tected in the o ptical (Buhrk e et allll988 | ). in the infrared (IR) 
( Reipurt h et alll2002l) . and in the radio ( Rodriguez & Reipurth, 
1 19961) . with a molecular cloud surrounding it. 

The HH 111 infrared source drives another highly colli- 
mated, long jet but the sourc e and the innermost 2 2" of the 
jet are very highly obscured ( Reipu rth et ail Il997l) . The ab- 
sorbed portions of the jet have been revealed in the infrared by 
Reip urth et all dl999l) and Reipur th et alJ (E oOO). Radio detec- 
tio ns have also been repo rted bv lRodrfguez*& Reipurth] dl994l) 
and[Reipurth et al. ( 1999). They identify a quadrupolar outflow, 
suggesting that the driving source is a close binary. 

3. Results 

3.1. Optically revealed, jet-driving T Tau stars 
3.1.1. Spectral interpretation 

Fig.dpresents EPIC PN spectra from DG Tau A, GV Tau, and 
DP Tau, respectively, together with spectral fits (see below). We 
have used the full exposure times except for a few intervals with 
high background radiation. The MOS spectra are similar but 
considerably fainter. They were used for spectral fits together 
with the PN spectra in the case of DG Tau A, while they were 
too faint to add useful information in the other cases. 

All three spectra are rather flat and show two shallow max- 
ima, one around 0.65-0.8 keV and the other at 1.5-3 keV, 
with an intervening trough at 1-1.5 keV. These spectra can- 
not be acceptably fitted with a combination of thermal coro- 
nal models subject to a single photo electric absorption compo- 
nent, as was already pointed out bv IGiidel et alJ ll2005h for the 
Chandra ACIS spectrum of DG Tau A. The softer maximum is 
formed by emission lines of O vm Lya at 18.97 A (0.65 keV), 
Fexvn at 17 A (0.73 keV) and at 15 A (0.83 keV), as well 
as Fexvm at 16 A (0.78 keV, blended with Ovm Ly/J) and 
at 14.2 A (0.87 keV). The Neix He-like triplet contributes at 
~ 13.5 A (0.92 keV, blended with Fxvn at 13.8 A or 0.9 keV). 
All these lines have maximum formation temperatures r max be- 
tween 3.2 MK and 6.9 MK (Ovm: 3.2 MK; Neix: 3.9 MK; 
Fexvn: 5.2 MK; Fexvm: 6.9 MK). The best matches of the 
subpeaks at 0.7 keV and at 0.80-0.85 keV in the DG Tau A 
spectrum are clearly the Ovm Lya and the two Fexvn lines. 
On the other hand, we note the absence of strong emission 
around 1 keV usually presen t in very active stars with hot 
coronae JTelleschi et allEo05l) . This latter emission would be 
formed by Ne x at 12.1 A (1.0 keV, T max = 5.9 MK), Fe xx at 
12.8 A (0.97 keV), and by various further lines of Fexx-xxiv 
(Tmax ~ 12-20 MK). For typical elemental compositions, 



these spectral properties clearly suggest emission by a rather 
cool plasma in the absence of prominent contributions from hot 
components. The rather high emission level down to 0.2 keV in 
DG Tau A and to a 0.3 - 0.4 keV in DP Tau further suggests a 
rather low absorbing column density. 

On the other hand, the spectra above » 1 .5 keV are shallow 
and can, despite the low count rates, be followed out to 10 keV 
for all three stars. This unambiguously indicates emission from 
very hot (> 10 MK) plasma which would normally form strong 
emission lines around 1 keV as well. The hard emission is dom- 
inated by bremsstrahlung which, in the absence of absorption, 
rises toward lower energies. But these low-energy features of 
hot plasma are suppressed by strong photoelectric absorption: 
All three hard spectral components reveal a turnover around 
2 keV. 

As suggested by iGiidel etafl J2005I) . these spectra can be 
successfully fitted with a model consisting of two thermal com- 
ponents, each subject to a separate photoelectric absorption 
component. We fitted the spectra as follows. For DG Tau A and 
GV Tau, we first considered counts within the soft maximum 
only (0.25-1.0 keV). The abundances adopted (and held fixed) 
reflect an "inverse First Ionization Potential" distribution often 
found in magnetically active and pre-main sequence stars (see 
IGiidel et alJ2006l for a discussion). Fit parameters were the ab- 
sorbing hydrogen column density A^h.s, a single temperature T s , 
and the emission measure EM S . The fits were perfect, with re- 
duced x 1 ~ 1 ■ Next, we added the hard portion of the observed 
spectrum up to 10 keV but held the fit obtained above fixed. 
We added a second spectral component with its own absorption 
component and fitted the corresponding variables A^.h, T^, and 
EMh. To optimize in particular the overlapping spectral region 
(1-1.5 keV), we continued fitting the A^h, T, and EM param- 
eters of both components simultaneously, but this did not sig- 
nificantly change the parameters. The DP Tau spectrum is too 
poor to allow for step- wise fitting. We therefore fitted the two 
components simultaneously from the outset. 

As far as abundances are concerned, the spectrum of DG 
Tau marginally suggested a somewhat higher Fe abundance 
(0.4 instead of the adopted 0.2 times th e solar photospheric 
value given by lAnders & Grevessell 19891) . but this value was 
poorly constrained, and the other parameters of interest here, 
kT and A^h, did not change in any significant way. We further 
experimented with variable abundances as follows. We kept the 
adopted abundance ratios fixed but varied the absolute abun- 
dance level. For DG Tau A, the best fit resulted in unaccept- 
ably low abundances of Fe (< 1% of the solar value), while the 
temperatures and absorption column densities again remained 
similar. 

The final results are reported in Table [3] The X-ray lu- 
minosities Lx,s and Lx,h refer to the 0.1-10 keV range based 
on integration of the best-fit model, an d assume a distance 
of 140 pc of the Ta urus association (e.g.. lLoinard et al1 f2005: 
iKenvon et alll994l) . 

Noting that the A^h value for the soft component of DG 
Tau A are lower than those anticipated from optical extinction, 
and A^h is only marginally in agreement with the lowest Ay re- 
ported for GV Tau A, we next adopted a fixed A^h value for the 
soft component (marked by "=" in Tabled}, using the standard 
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DG Tau A GV Tau A DP Tail 




^0.2 0.5 1 2 5 10 ^0.2 0.5 I 2 5 10 ^0.2 0.5 1 2 5 10 



channel energy (keV) channel energy (keV) channel energy (keV) 

Fig. 1. Average spectra of DG Tau A (left), GV Tau (middle), and DP Tau (right). Also shown are the fits to the spectra (black 
histograms) and separately the soft (red) and the hard (blue) spectral components. 

Table 3. X-ray properties of the TAX sources 



Soft component Average hard component 



Star 


N" 








1 h 




* 2 /dof 




(10 22 cm- 2 ) 


(10 6 K) 


(10 29 ergs-') 


(10 22 cm" 2 ) 


(10 6 K) 


(10 29 ergs" 1 ) 




DG Tau A c 


0.11 (0.08-0.14) 


3.7 (3.2-4.6) 


0.96 


1.80(1.21-2.40) 


69 (44-188) 


5.1 


89.9/93 




= 0.45 


2.2 (2.05-2.35) 


9.6 


1.71 (1.14-2.30) 


73 (45-220) 


5.0 


100.2/94 


GV Tau A 


0.12 (0.01-0.79) 


5.8 (2.1-9.1) 


0.54 


4.12(2.07-6.47) 


80 (34-...) 


10.2 


20.3/19 




= 1.12 


0.94(0.94-1.31) 


4170 


3.48 (1.69-5.48) 


95 (42-...) 


9.6 


24.6/20 


DP Tau 


0.00(0.00-1.05) 


3.2 (0.0-7.7) 


0.04 


3.78 (1.38-14.0) 


61 (10-...) 


1.1 


8.9/7 




= 0.29 


2.3 (0.94-3.8) 


1.1 


3.13(0.78-13.8) 


104 (10-...) 


0.92 


9.1/8 



a 90% confidence intervals in parentheses (ellipses indicating unconstrained parameter) 

b Modeled for the 0. 1-10 keV energy interval 

c Fits applied to combined PN, MOS 1 , and MOS2 data 



conversion formula, Nh ~ 2 x 10 21 Ay cm" 2 , valid for stan- 
dard gas-to-dust ratios (see IVuong et alJl2003l and references 
therein). These fit results are also reported in Table Fixing 
TVh at the expected values provides significantly worse fits and 
lowers the temperature of the cooler component, at the same 
time rising its Lx- The value of the latter is unacceptable for 
GV Tau because it approaches L„, while Lx is usually bound by 
Lx/L, ^ 10~ 3 (Gudel 2004 and references therein). We believe 
that the origin of this excess is purely numerical: increased TVh 
allows more and cooler plasma to be present without significant 
spectral contribution, because X-rays from cooler plasma are 
softer and therefore more strongly absorbed. The effect contin- 
ues as one increases Nh- This lends most credibility to the best 
fit without enforced TVh, and this best fit indeed results in only 
modest suppression even at low energies. 

In summary, the spectra of DG Tau A, GV Tau, and DP 
Tau are thus composed of two components, a cool component 
subject to very low absorption and a hot component subject to 
photoelectric absorption about one order of magnitude higher. 
The cool component shows temperatures atyp ical for T Tau 
stars ( Skinner et al., 2003; Telleschi et al., 2006), ranging from 
3-6 MK (best fit), while the hot component reveals extremely 



high temperatures (60-100 MK). If Nh is fitted to the spectrum, 
it tends to be lower than expected from the optical extinction. 

The spectrum of CW Tau resembles the spectrum of DP 
Tau but is too weak to derive definitive results. This source 
is also affected by the wings of the nearby and very bright 
X-ray source of V773 Tau (at a distance of 1.5'). A single-T 
fit to the harder photons (above 1.2 keV) requires an exces- 
sively absorbed (Nh ~ 10 23 crrr 2 ), very strong cool compo- 
nent with T « 4 MK, but there is a marginal excess of counts 
below 1 MK that requires a further, low-absorption compo- 
nent. Better-qua lity data are needed. We no te that the star's 
Ay = 2.29 mag ( K en von & Hartmannl 1 19951) suggests a much 
smaller Nh ~ 5 x 10 21 cirr 2 . 

Results for the three further jet/outflow driving sources 
mentioned above, UZ Tau, D P Tau, and XZ Tau, are sum- 
marized in iGiidel etafl d2006l) . UZ Tau (a triple system) re- 
veals a relatively soft spectrum (emission-measure weighted 
average temperature of 5-7 MK). On the contrary, the very 
faint spectrum of DD Tau is extremely hard, reminiscent of the 
harder component of the sources discussed above, with T- dv = 
38-42 MK and excess absorpti on (7V H = T2.3-2 91xl 21 cirr 2 , 
whereas A v = 0.39 mag after IWhite & Ghe dEoOll or Aj = 
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0.05 mag after lBriceno et alJl2002h . A soft excess is not found 
for this source. Adopting a jet bulk vel ocity of 300 km s -1 and 
a radial velocity of v ra( j at 75 km s _1 dHirth et all Il997l) . we 
infer an inclination angle of »75 deg. The strong photoelectric 
absorption may be due to disk material, while the low values of 
Ay and Aj may be affected by scattered light. No anomaly was 
found for XZ Tau. 

3.1 .2. Light curves and time-resolved spectroscopy 

We present the light curves in Fig. |2] where we used counts 
below at 1 keV and above « 1.6-1.8 keV (depending on 
the spectral shape) for the soft and the hard component, re- 
spectively. All light curves show clear signatures of flaring. 
Surprisingly, however, flaring is revealed only in the hard com- 
ponent while the soft component remains at a low level. We 
tested the background-subtracted light curves for constancy, 
applying a x 1 t est - To do so, we rebinned the light curve of 
DG Tau A such that the bins of the soft curve contained, on 
average, 8 or 11 counts. This required bins of 1500-2000 sec- 
onds, which is sufficient to still recognize typical flares. With 
the same binning, the hard curve contained, on average, 6 or 
8 cts but showed the strong flare well resolved. For GV Tau, 
we binned to 2000 and 3000 s, which on average resulted in 
only 2.3 and 5 cts per bin, respectively. The hard curve shows 
a 2.7 times higher count rate. A test for constancy resulted in a 
reduced y 1 that was smaller than unity for the soft light curves 
of DG Tau A and GV Tau, whereas the variability in the hard 
curve is highly significant with^ d » 10- 17. DP Tau is too 
faint for a meaningful^- 2 test of this kind. 

To further clarify the spectral behavior in time, we extracted 
EPIC spectra during various episodes for DG Tau A and GV 
Tau (low level, flare rise, and flare decay for DG Tau A; low 
level, first flare episode and second flare episode for GV Tau; 
see Fig. |2] for a graphical definition of the corresponding time 
intervals). The corresponding plots are shown in Fig. [3] The 
hard component varies both in flux and in spectral shape: the 
shallower slope during the flares indicates higher temperatures, 
while the soft component remains at the same level. This is 
borne out by spectral fits performed for the same intervals. The 
results are reported in Table|4] Because the spectra are too poor 
to fit Nh reliably, this parameter was frozen at the value found 
for each integrated spectrum. The temperature of the soft com- 
ponent remains constant within the errors, while the tempera- 
ture of the hard component rises during the flare episodes. Note 
also that the X-ray luminosities of the soft component are very 
similar (the spectrum of the first flare episode of GV Tau is too 
faint to produce reliable results). 

3.1 .3. Comparison with Chandra spectra 

Both DG Tau A and GV Tau have been observed by Chandra 

I II I 

ACIS as well (see Giidel et al. 2006 for an observation log). 
The sources are rather faint, and no flaring was detected during 
the observations. We therefore present only the integrated spec- 
tra and compare them with the integrated XMM-Newton spec- 
tra. Figure@]shows the observed Chandra spectra together with 
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Fig. 2. Background-subtracted light curves of DG Tau A (top), 
GV Tau (middle), and DP Tau (bottom), binned to 500 s, 500 s, 
and 1000 s, respectively. Each figure contains three panels, 
showing, from top to bottom: Soft photons (0.4-1.0 keV for 
DG Tau A and GV Tau, and 0.35-1 keV for DP Tau); hard 
photons (1.6-7.3 keV for DG Tau A, 1.7-7.3 keV for GV Tau, 
1.8-10 keV for DP Tau); full band (0.4-7.3 keV for DG Tau A 
and GV Tau, and 0.35-10 keV for DP Tau). For DG Tau A and 
GV Tau, intervals separately studied are marked by red vertical 
lines. Only the time intervals observed jointly by all three EPIC 
cameras are shown. 
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Soft component Average hard component 



Star 








w 


x h 


L X.h 


* 2 /dof 




(10 22 CITT 2 ) 


(10 6 K) 


(10 29 erg s" 1 ) 


(10 22 cm" 2 ) 


(10 6 K) 


(10 29 erg s _1 ) 




DG Tau A low c 


=:0.11 


3.7 (3.3-4.6) 


0.91 


= :1.80 


23 (14-42) 


2.1 


24.5/33 


DG Tau A rise 1 


=:0.11 


3.9 (3.1-5.6) 


0.89 


=:1.80 


396 (144-...) 


12.4 


27.9/31 


DG Tau A decay c 


= :0.11 


3.9 (3.3-4.9) 


1.02 


= :1.80 


74 (44-144) 


6.7 


15.7/24 


GV Tau low 


=:0.12 


4.7 (3.0-9.7) 


0.54 


:=4.12 


17.5 (4.3-63) 


7.3 


5.4/5 


GV Tau flare 1 


=:0.12 


2.1(...) 


0.16 


= :4.12 


77 (24-...) 


16.3 


3.7/2 


GV Tau flare 2 


= :0.12 


4.5 (...-12.6) 


0.45 


= :4.12 


119(52-...) 


19.9 


6.7/7 



a 90% confidence intervals in parentheses (ellipses indicating unconstrained parameter) 

b Modeled for the 0. 1-10 keV energy interval 

Fits applied to combined PN, MOS 1 , and MOS2 data 
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Fig. 3. Upper panel: Spectra and fits of DG Tau A during low-level episodes (left), during the rise to the large flare (middle) and 
during the decay phase (right). Lower panel: Spectra and fits of GV Tau during low-level episodes (left), during the first flare 
episode (middle) and during the second flare episode (right). 



the spectral model of the soft component of the XMM-Newton 
observations, folded with the appropriate ACIS response matri- 
ces. We note that the Chandra and XMM-Newton observations 
were separated by 7 months for DG Tau A and by 8 months for 
GV Tau. 

Whereas the flux of the hard components varied greatly 
during the 30-40 ks XMM-Newton observations, the flux of the 
soft components changed little in the 7-8 months between the 
Chandra and XMM-Newton observations, although some de- 



tails in the spectrum do differ. The Chandra spectrum of DG 
Tau A shows excess flux at 0.8-0.9 keV, a discrepancy that did 
not disappear if the the soft portions (<1 keV) of both spectra 
were fitted simultaneously, even if we changed the abundances. 
Fitting the Chandra spectrum find s a best-fit Fe abundance of 
» 0.2 (Ander s & GrevesseL 1 1989). the same that we adopted 
for the XMM-Newton model. The soft-component 7V H is again 
low(« 0.6 xlO 21 cm" 2 ) and marginally compatible with XMM- 
Newton, although the best-fit temperature is somewhat higher 
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DG Tau A, Chandra 
^_ Average spectrum 
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Fig. 4. Chandra spectrum of DG Tau A (from ACIS-S) and GV 
Tau (from ACIS-I). The histogram corresponds to the soft com- 
ponent modeled from XMM-Newton data but folded through 
the appropriate ACIS responses. 

(4.4 MK). No significant deviation between the soft Chandra 
and XMM-Newton spectra is found for GV Tau. 

3.2. Embedded, jet-driving protostars 

The embedded, jet-driving protostars seen at high inclination 
expectedly reveal rather faint X-ray fluxes owing to strong 
photoelectric absorption. Table [5] summarizes X-ray detection 
properties and results from spectral fits (see below). 

Three of the four target sources (DG Tau B, FS Tau B, and 
HH 34 IRS) were clearly detected very close to the expected 
positions based on wavdetect, while only 4 counts were reg- 
istered for HH 111 (Fig. 0. The close positional coincidence 
with the radio s ource (0.15" distance from th e HH 111 radio 
position given by Rodriguez & Reipurth 1994), the strong clus- 
tering of the counts within a radius of 1-2 pixels (Fig.[5J, the 
energy of the four counts (3.3-6.4 keV, very similar to the dis- 
tributions found for the other absorbed sources), and the rather 
low background level make it probable that these photons orig- 
inate from the HH 111 jet-driving protostar although the signif- 
icance is too low to prove this identification. The background 




■ 30" 



hh ml 




Fig. 5. Chandra ACIS-S images of jet-driving sources. Pixel 
size is 0.49"; all images were exposed for a 30 ks and show 
counts in the 0.5-7 keV range except for HH 111, where only 
counts in the 3-7 keV range are shown to minimize back- 
ground. From top to bottom: DG Tau A and B (separation 
* 50"); FS Tau A and B (separation * 20"); HH34 IRS; 
HH1 1 1 IRS (tentative detection). 

level was such that 0.1 counts would be expected in the 3- 

7 keV range in a circle with a radius of 3" around the source, 
and 0.6 cts in the 0.3-10 keV range. 

To avoid systematic effects in the spectral fits due to coarse 
binning of the few counts, we used the C statistic in XSPEC 
in combination with unbinned data in the energy range of 0.5- 

8 keV. The metallicity was fixed at 0.2 times solar photospheric 
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Table 5. Protostellar, jet-driving X-ray sources 



Parameter 


DG Tau B 


Ref 


FS Tau B 


Ref 


HH 34 IRS 


Ref 


HH 1 1 1 IRS 


Ref 


Useful spectral range (keV) 


2-6.5 




2-7 




3-7 




3-6.5 




Detected X-ray counts 


9 




32 




15 




4 




Position: 


















Expected R.A. (h m s) 


04 27 02.55 


1 


04 22 00.70 


2 


05 35 29.84 


2 


05 51 46.31 


3 


Expected 6 (deg ' ") 


26 05 30.90 


1 


26 57 32.50 


2 


-06 26 58.40 


2 


02 48 29.72 


3 


Observed R.A. (h m s) 


04 27 02.58 




04 22 00.71 




05 35 29.84 




05 51 46.30 




Observed 6 (deg ' ") 


26 05 30.91 




26 57 32.17 




-06 26 58.24 




02 48 29.70 




Offset (") 


0.40 




0.36 




0.16 




0.15 




X-rays: 


















Expected N 1 ^ 


1.4 


6 


> 1.6, £ 2,4.6 


1,6,4 


1 


5 


> 6 


7 


Measured 


51 (6-72) 




26 (10-37) 




28 (14-54) 









" References: 1 = Padgett et al. (1999), 2 = 2MASS, Cutri et al. (2003), 3 = Rodriguez & Reipurth (1994), 4 = Krist et al. (1998), 
5 = Reipurth et al. (1986) . 6 = Mundt et al. (1987) . 7 = Reipurth (1989) 
b Expected N a from reported optical extinction: Nn ~ 2 x 10 21 A V cirr 2 
Njj is given in units of 10 22 cirr 2 , 90% confidence ranges are in parentheses 



DG Tau B observed (black), unabsorbed (red) FS Tau B, observed (black), unabsorbed (red) HH34 IRS, observed (black), unabsorbed (red) 



o 




channel energy (keV) channel energy (keV) channel energy (keV) 



Fig. 6. Binned ACIS-S CCD spectra (red/gray, with error bars) of (from left to right) DG Tau B, FS Tau B, and HH34 IRS, 
compared with the modeled, response-folded intrinsic spectra if the photoelectric absorption is removed (black solid histograms). 



values ( Anders & Grevesse Il989i) . coincident with the adopted 
Fe abundance for the other spectral fits (we note that the only 
emission lines significantly contributing to the absorbed pro- 
tostellar sources are those in the Fe complex at 6.7 keV; even 
adopting solar abundances did not change our results signif- 
icantly). For illustration purposes, we show in Fig. [6] binned 
spectra together with the best-fit unabsorbed models. These 
figures illustrate that most of the flux is suppressed by photo- 
electric absorption. The best-fit results and the acceptable 90% 
ranges for A^h are summarized in Table |5j other parameters 
were not sufficiently constrained for further discussion. 



4. Discussion 

The TAX spectra from the optically revealed jet-driving T Tau 
stars presented in this paper clearly require two spatially dis- 
tinct sources because both their time behavior and their ab- 
sorption along the line of sight are different. We now discuss 
possible models. 



4.1. CTTS: The hard component 

The hard components of the spectra reveal very high tempera- 
tures. Such tempera t ures a r e commonly found i n other T Tau 
stars ( Skinne r et all l2003t iTelleschi e t al.. 2006). The flares 
seen in the light curves clearly argue in favor of a magnetic ori- 
gin. For DG Tau A, the U band light curve extracted from the 
simultaneous observations with the OM peaked before the X- 
rays, in analogy to solar flares (Fig. [7] upper panel). We there- 
fore interpret this component as being due to a magnetically 
confined corona above the surface of the star. The time reso- 
lution of the OM data is too coarse to reveal a clear temporal 
relation for GV Tau (Fig. [7] lower panel). 

The excessive absorption of this coronal component re- 
quires a large amount of cool gas along the line of sight to 
the observer. A cool wind or a molecular outflow are potential 
absorbers. However, in that case the soft component must be 
located at very large distances from the star in order to escape 
from strong absorption. Also, dust admixtures would lead to 
correspondingly strong optical extinction of the star, which is 
not observed. 

We propose a scenario related to accretion streams. The 
mass accretion rates of our stars are among the largest found 
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DG Tau A 




Haro 6-10 AB 



3242.90 




3243.10 3243.20 

JD - 2450000.0 



Fig. 7. Correlation between X-ray emission and U band emis- 
sion. The U band data are shown by the large crosses above 
the X-ray light curves and refer to the y axis labels on the right 
side. Top: DG Tau A; bottom: GV Tau. 

in any TTS in the TMC (Table |5Jl. Much of the gas accreting 
from the circumstellar disk is only weakly ionized. This gas 
will stream along the magnetic field lines that form the corona 
and thus enshroud X-ray bright coronal loops so that photoelec- 
tric absorption attenuates the soft X-rays from the underlying 
coronal plasma. 

Why, then, is the optical extinction of the stars relatively 
small? The measured Nh columns would imply Ay » 10 — 
20 mag for each star while measured optical extinctions are 
only a few magnitudes (Table |2jl. Such deviations are possi- 
ble if the X-rays are propagating through an anomalous mix- 
ture of gas and dust in which dust is depleted. Dust destruc- 
tion and evaporati on occurs at temperatures around T sn h ~ 
1600 - 2000 K to'Alessio et all Il998l IWhitnev et all l2004 . 
and the dust sublimation radius from the center of the star can 
be estimated from the empirical formula, 



(1) 



"sub _ I 1 sub I 

R* \ T e f[ I 

JWhitnev et alJ2004l where T m \, = 1600 K has been assumed). 
Using T e ff from Table[2] we find radial distances of 7 - 10 R*. 
As a consequence, there should be no strong extinction by dust 
from the accretion streams while photoelectric absorption re- 
mains in effect. In this picture, visual and near-infrared extinc- 



tions are due to the larger-scale environment of the star, while 
photoelectric absorption is due to both the large-scale environ- 
ment (for the soft components) and, more strongly, the imme- 
diate surroundings of the star (for the hard component). Our 
observations may thus give indirect evidence for dust sublima- 
tion in the inner circumstellar accretion disks of T Tau stars on 
the one hand and for strong gaseous accretion streams in the 
innermost region or onto the star on the other hand. 

4.2. CTTS: The soft component 

The soft component dominates the spectrum below 1 keV. The 
spectral shape in this region shows no indication of an equally 
prominent plasma component at temperatures ^10 MK. This 
component, considered in isolation, is unusual for T Tau stars 
that normally reveal spectra dominated by hot temperatures 
jTelleschi et allE()06l) . The exception extensively discussed in 
the recent literature is TW Hya, the X-ray spectrum of which is 
strongly dominated by a soft spectral component of this kind. 
TW Hya is, however, a rather evolved classical T Tau star with 
an age of » 1 Myr and a low accretion rate of 4x 1 0~ 10 M Q yr~' 
(Muz erolle et allboOCfc . lKastner et all 12002) proposed that the 
low temperatures and the high electron densities inferred from 
the O VII He-like triplet are indicative of X-rays produced in 
accreti on shocks rather than coronal X-ravs. lstelzer & Schmittl 
(2004) argued that the anomalous abundances found in the TW 
Hya X-ray spectrum are further support for this picture. 

In the case of the jet-driving sources, this explanation faces 
serious difficulties. The accretion rates measured for these stars 
are among the highest known in TMC. If X-rays are produced 
near the footpoints of accretion funnels, then we would expect 
that the same gas streams absorbing the hot coronal emission 
also strongly absorb the soft component. 

On the contrary, the measured absorption column densities 
tend to be smaller than expected from the optical extinction 
of the star, at least for DG Tau A. Either, the soft X-rays es- 
cape through a region of depleted gas in between the accretion 
streams, or they are formed significantly outside the immediate 
environment of the star. Consulting Table|2] the distinguishing 
property of the TAX sources are their jets. The strong differ- 
ences in Lx of the cool component (ratio of 1 :0. 56:0. 04 for DG 
Tau A : GV Tau : DP Tau) mirror in two parameters: the equiv- 
alent width of [O I] (1:0.23:0.05), and the mass outflow rate 
(1:0.3:0.06). Although no linear relation is to be expected, the 
trend is rather suggestive for outflow-related X-ray activity. 

A faint soft component was also found at a distance of 2- 
5" from DG Tau A, co -spatial with bow-shock structures in the 
jet dGudel et alJ2005l see also Fig.|5j. Because the knots at this 
dist ance show X-ra y emission while those further away do not, 
iGtidel et"all (l2O05) proposed that the more luminous soft com- 
ponent forms in shocks even closer to the star, presumably at 
the base of the jets but clearly outside the magnetospheric ac- 
cretion zone around the star. This would simply require that the 
jet portions closer to the star contain more emission measure 
heated to X-ray temperatures although the temperatures would 
be the same. This scenario is supported by the close similar- 
ity between the soft component of the X-ray spectrum obtained 
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from the region within one arcsecond around the star, and the 
extended X-ray source coincident with the jet. 

Further support for this hypothesis comes from optical ex- 
tinction measu rements of the Herbig -Haro structures in the jet 
of DG Tau A. ICohen & Fuller! (Il985h found an Ay of 0.39 mag 
while the stellar A\ has been reported between 1 .4 mag and 
3.3 mag (Table |3. Using standard interstellar gas-to-dust ra- 
tios, we would expect a hydrogen column density of A^h ~ 
8 x 10 20 cm~ 2 toward the jet, comparable with the X-ray de- 
termined value of (1.1 + 0.3) x 10 21 cirT 2 . 

The precise ca use for the heating of such jet/related X-ray 
sources is unclear. iBallv et al] J2003I) suggested, for an X-ray 
source in the jet but spatially resolved from the location of the 
fully absorbed protostar L1551 IRS-5, a number of scenarios: 

i) Shock-heating at the working surface of a bow shock col- 
liding with the ambient medium, in analogy to shock heating in 
the more distant HH objects; ii) Thompson scattering of proto- 
stellar X-rays into the line of sight by a cool, dense medium 
perhaps in the outer, expanding regions of the circumstellar 
disk or the accreting envelope; iii) shock heating at an obsta- 
cle immersed in the jet flow, where either magnetic fields or a 
dense, ambient medium redirects the initially expanding wind 
and collimates it to a jet; iv) wind-disk or wind-wind collisions 
in a binary system. We are not in a position to discriminate 
between such models, except to note that the final option is un- 
likely to apply to our sources, as discussed in the next section. 

We speculate that magnetic fields may further heat the 
jet gas. The jet of DG Tau A has been found to rotate 
around its axis. This motion derives from the Kepler mo- 
tion of the circumstellar disk from where the jet acceler- 
ates. The motion is differential in radius, an d the mass flow 
along the field generates a trailing spiral \ Anderson et all 
2003J. In ionized gas, such field configurations induce cur- 
rents that may dissipate and heat the gas if its resistivity is 
large enough. Further irregularities can be introduced by pulsed 
outfl ows, precession, and the shock with t he ambient medium 
JCera ueira & de Gouveia Dal Pintol l2004). Magnetic "tangen- 
tial discontinuities" may lead to reco nnection, thus liberating 
additional (magnetic) energy llParkeii fl983>. Direct magnetic 
field measurements are yet to be made, although radio obser- 
vations sugge st the presence of weak magnetic fields in jets or 
outflows (e.g., iRav et alJl997l) . 

We next compare in Fig. [3] three CCD spectra from 
Chandra (ACIS-S): The spectrum of DG Tau A (black, with 
soft excess), the spectrum of FS Tau A (thick, blue), and the 
spectrum of the protostar FS Tau B. 3 FS Tau A is subject 
to a similar absorption to the h ot component in D G Tau A 
(N H = [0.9 - 1.4] x 10 22 cm' 2 . Idiidel e"tal] E)06). shows a 
similar (unabsorbed) Lx ~ 3.4 x 10 29 erg s , and also reveals 
an extremely hot electron temperature of 43 MK (higher tem- 
perature of a two-component fit; Gud el et alfZOOfih . FS Tau is, 
however, not known to drive strong j ets, although a poo rly col- 
limated wind seems to be present IHirth e t al.. 1997). While 

3 we discuss data only from Chandra for this comparison because 
the strongly absorbed FS Tau B is located in the wings of the XMM- 
Newton image of FS Tau A, and FS Tau A itself was showing a gradual 
decay from a strong flare during the XMM-Newton observation 
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Fig. 8. Comparison of the three Chandra spectra from DG Tau 
A (black, with soft component), FS Tau A (blue, thick), and FS 
Tau B (red). 



the two spectra are very similar above 1 .5 keV, the absence of 
a soft component in FS Tau A is striking. The spectrum of the 
jet-driving FS Tau B does not show any soft emission because 
such emission from the star and the jet within an arcsecond or 
so would be absorbed. The hard component above 5 keV, how- 
ever, is again similar in flux to the hard spectral component of 
the other two stars, suggesting a similar coronal component. 

Optical jets are detected in predominantly low-excitation 
lines such as forbidden lines of [O I], [N III , and [Si III forme d 
at no more than a few thousand degrees llHirth et all fl99l . 
On the other hand, the X-ray emission de tected in direct imag- 
ing around DG Tau A (Giidel et al., 2005) and the soft spectral 
components suggested here to be due to jets as well require 
plasma with electron temperatures above 1 MK. The challenge 
is to reconcile the observational evidence from two widely dis- 
parate temperature ranges when the presence of an intermediate 
temperature range ( 1 4 - 1 6 K) is apparently less evident in the 
observations. What is the status with respect to intermediate- 
excitation lines? 

First, detectability of forbidden lines requires favorable 
densities for their formation, and a sufficient emission measure. 
Modeling of [O I], [N I II, and [Si III lin es indicates that they 
form in distinct regions (H am an rl 11994 : [O I] forms in a gas 
in the electron-temperature range of 9,000-14,000 K and high 
densities (n e « 5 x 10 5 - 10 7 cm -3 ), while [Si II] forms above 
13,000 K and in lower densities (n e * 10 3 - 7 x 10 4 cm" 3 ). 
At high velocities only, [N II] and [O II] lines may be formed, 
requiring r e <: 15,000 K and n e 10 s cm 4 . These results 
therefore indicate that n e decreases away from the star, while 
T e increases, at least as far as high-velocity gas is concerned. 

Conditions in the jets may in general simply not be favor- 
able to the production of detectable amounts of flux in higher- 
excitation lines, e.g., [O III]. On the other hand, there are now 
several rep orts of the presence of hot winds and outflows from 
T Tau stars. Beristain et al. (2001) presented an extensive study 
of He I and He II emission line profiles. The important diagnos- 
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tic power of these lines lies in their high excitation potentials 
(20-50 eV), requiring gas tempera tures of (2,5 - 9) X 10 4 K in 
the case of collisional excitation jBeristain et all 1200 ll) . The 
authors found that blueshifted absorption features in broad 
He I line profiles and maximum blue-wing velocities exceed- 
ing 200 km s _1 indicate the presence of hot, coronal winds in 
about half of their sample, one of the clearest examples being 
DG Tau A with a maximum blue- wing velocity of 600 km s _1 . 
These winds would be launched in the polar region of the stars, 
complementary to cooler winds that are probably accelerated 
around the disk/magnetosphere boundary. 

Work on He line s was extended t o the H e I /110830 line by 
ITakami et all J2002h . lEdwards etalJ J2003I) and lEdwards et all 
(2006). This feature is again a high-excitation line (20 eV, re- 
quiring T * 10 4 K for collisional ionization) and has an un- 
precedented sensitivity to hot, inner winds. Significant absorp- 
tion below the continuum was found in 71% of the CTTS, 
but in none of the WTTS, suggesting that the feature is due 
to accretion-driven winds, originating in the immediate vicin- 
ity of the star. Specifically in the case of DG Tau A, ab- 
sorption and emission components were found, identified with 
an in ner wind and the jet shocks, respectively ( ITakami et all 
120021) . Further, [O III] ,14959 and A5001 emission is indeed 
also foun d in one of the ou ter working surfaces of the DG 
Tau jet JCohen & Fuller!! 19851) and in other shock-r elated fea- 
tures in Herbig-Haro objects Jlvlatt & Bohmll2003l and refer- 
ences therein), indicating that singly ionized O is further ion- 
ized beyond its ionization potential of 35 eV. Shock model- 
ing suggests post-shoc k electron tempera tures in the range of 
(6 - 7) x 10 4 K llCohen &Fullei!ll985l) . and higher temper- 
atures at specific locations in Herbig-Haro shock regions must 
exist to explain recently detected X-ray emission (see summary 
in Sect. [Hi. 

iGomez & Verdu^ol J200lh have analyzed high-resolution 
profiles of semiforbidden ultraviolet lines of C III] A 1908 and 
Si III] ,11892 in two CTTS. They concluded that the lines are 
not formed in accretion shocks but in a hot wind outside 2R* 
and not farther from the star than 38/?* (for RY Tau), suggesting 
a bow-shaped shock source at the base of the jet, with densities 
> 10 5 cm~ 3 (i.e., at least 10 2 - 10 3 times denser than Herbig- 
Haro objects) and temperatures in the range of 5 x 10 4 - 10 5 K. 

Still higher outflow temperatures have been suggested by 
far-ultraviolet observations. The CTTSs TW Hya and T Tau 
both show P Cygni profiles, line asymmetries, and absorption 
in C III ,1977 and O VI ,11032 lines that are indicative of a fast 
(400 km s" 1 ) and hot (3 x 10 5 K) accelerating outflow close to 
the star dDupree et alll2005l) . 

In summary, although there is no coherent model of winds 
and jets at all temperature and density levels, it appears clear 
that CTTS show evidence for jet- or wind-related gas at tem- 
peratures from < 10 4 K to several times 10 5 K. It is then natural 
to think of the X-ray emission as a continuation toward higher 
temperatures. It appears that at various locations in the jet, in 
particular in shock regions, conditions may be favorable for the 
production of X-rays, while densities and volumes seem to be 
unfavorable for appreciable line emission at lower ionization 
stages such as [O III]. A speculative possibility would be that 
rapidly heated X-ray sources in the jets freely expand and cool 



Table 6. Coordinates of GV Tau A, B 



Object 


Observation 


RA(J2000.0) <! 
h m s 


5(J2000.0) fl 
deg ' " 


GV Tau A 
GV Tau B b 


2MASS 
2MASS 


4 29 23.734 
4 29 23.725 


+24 33 00.277 
+24 33 01.577 


GV Tau A c 
GV Tau B c 


VLA 
VLA 


4 29 23.733 
4 29 23.724 


+24 33 00.125 
+24 33 01.425 


softX 
hardX 
total X 


Chandra 
Chandra 
Chandra 


4 29 23.744(8) 
4 29 23.749(7) 
4 29 23.737(4) 


+24 33 00.54(11) 
+24 33 00.66(10) 
+24 33 00.53(5) 



a Errors for Chandra sources in units of last digits; formal wavde- 
tect fit errors 

b Appl ying radio offset to GV Tau A 2MASS source 
c From R eipurth et al J l200l) 



in such a way that their densities and emission measures are too 
low for easy detection in intermediate-excitation lines once the 
temperature has dropped to appropriate values. The dominant 
[O I] emission i s not necessarily formed at the same locations 
ilHamann![l994T) . 

4.3. Binarity? 

TAX spectra may originate from a binary with components that 
are located behind largely different gas columns, each thus con- 
tributing one of the spectral components. This model is un- 
likely for the following reasons: i) The less absorbed com- 
panion would in all cases reveal a uniquely soft, non-flaring 
X-ray component only. None of the bright TTS in the XEST 
survey, single or mul t iple, revealed such an X-ray spectrum 
llTelleschi et all 120061 Idiidel etall 120061) . ii) Except for GV 
Tau, no companions have be en found to the TAX sources de- 
spite detailed searches (e.g., Lei nert et al.lll99ll) . iii) GV Tau 
= Haro 6-10 is indeed a binary, the more absorbed compo- 
nent being an embedded protostar behind a large gas column 
(see lReipurth et al.ll2004l for a radio image). A high-resolution 
archival Chandra image (ObsID = 4498), however, reveals that 
the soft and hard photons originate from the same location, 
and this l ocation agrees, within the error ranges for Chandra, 
the VLA ( Reipurth etall \2Q04 and 2MASS, with the less ab- 
sorbed component Haro 6-10A (Table|6| Fig.|5)- We conclude 
that the embedded Haro 6-1 0B does not contribute to the de- 
tected X-ray emission. 

We also note that the prominent jet from the GV Tau system 
is associated with the optically revealed star G V Tau A r ather 
than the deeply embedded infrared source dReipurth et all 
120041) . 

4.4. The protostars 

Although the ACIS protostellar sources are very faint, their 
photon energy distribution is revealing: we detect at most one 
count below * 2 keV in each source, which is normally the 
dominant spectral range for a coronal source. 



Manuel Giidel et al.: X-rays from jet-driving protostars and T Tauri stars 



15 



■ 



Fig. 9. Chandra images of GV Tau. Left: Counts in the range 0.5-5 keV. The lower cross marks the 2MASS posi tion of GV Tau 
A, th e upper cross the GV Tau B protostar, applying the same offset to GV Tau A as measured in radio data ( Reipurt h et all 
2004). Middle: Similar, for the 0.5-1 keV range. Right: Similar, for the 2.2-7 keV range. Pixel size is 07492. 



The emission cannot be radiation from shocks in the jets or accretion shock close to the surfac e of the star as propo sed for 



Herbig-Haro objects. Such sources are too soft and will there- 
fore be absorbed tBallv et all2003l) . The hard sources required 
here are reminiscent of coronal emission of very active stars. 
Useful information on the hydrogen absorption column den- 
sity, A^h, can therefore be obtained/or the circumstellar disk. 

Although the 90% confidence ranges in Nu, kT, and Lx are 
large (Table [5jl, the values of A^h are consistently found in the 
range of » (5 - 50) x 10 22 cirT 2 if we require "reasonable" 
values for kT and L x (e.g., typical saturation values given the 
estimated L*, and T in the range of 6 - 60 MK). We thus expect 
Ay ss 30 - 30 mag, mu ch higher than reports in the literature 
(A v * 5 -25:lReipurth et alJl986lfReipurtrJl989HMundt etai] 
Il987l iKrist et aljll998l P991: Either, the A v values are erro- 
neous due to IR emission from the stellar surrounding, or the 
gas-to-dust (G/D) ratios are significantly enhanced in these cir- 
cumstellar disks. 



5. Summary and conclusions 

Jets are assumed to be launched within a few AU from the 
central star and have been resolved down to 0.1" (14 AU) 
from t he latter at th e dist ance of the Taurus Molecular 
Cloud ( iBacciotti et all 120001) . These regions remain unre- 
solved from the central star in present-day X-ray observations. 
Nevertheless, several jet-driving T Tau stars reveal peculiar X- 
ray spectra that have been ta ken as evidence for shock for- 
mation at the base of the jets JGudel et alJl2005t and our dis- 
cussion above). In all cases, a very weakly attenuated (A^h ~ 
10 21 cm -2 ), luminous but very soft spectral component indi- 
cates electron temperatures of no more than a few MK. While 
this component produces photons mostly between 0.5-1 keV, 
a strongly absorbed (Nh > 10 22 cm -2 ), hard spectral compo- 
nent dominates above 2 keV, resulting in a shallow spectral 
shape with two peaks. The hard component originates from hot 
plasma of several tens of MK. In each case, it reveals flares 
while no significant variability is detected from the soft compo- 
nent whose flux is approximately constant, even on time scales 
of 6-8 months (Sect. l3~T3t . 

We attribute the hard component to a magnetospheric 
corona, as observed on most T Tauri stars. Although these 
T Tauri stars are unusually strong accretors (with M acc up to 
tss 10~ 6 M yr -1 ) we do not attribute the soft component to an 



accretion snock close to tne surface of tne star as proposed tor 
the so far unique CTTS TW Hya jKastner et all 12002b . Such 
emission would be subject to the same high absorption experi- 
enced by the hard component, while the observed soft compo- 
nent is little absorbed. Instead, we propose that the soft compo- 
nent originates further away from the star where the first shocks 
form at the base of the jets. Indeed, the X-ray luminosity of the 
soft component in these three stars appears to scale with the 
mass outflow rate and with the equivalent width of [O I] from 
the jet. 

It is interesting to note that the visual extinction of the T 
Tau stars showing Two- Absorber X-ray (TAX) spectra is mod- 
est and does not agree with the photoelectric absorption of the 
hard component if standard gas-to-dust mass ratios compatible 
with the interstellar medium are assumed. We suggest that the 
hard emission is absorbed by accreting material an that this ac- 
creting material is dust-depleted. This is supported by the dust 
sublimation radius of our sources being several stellar radii. 
Our observations thus provide both evidence for dust sublima- 
tion in the innermost part of the accretion disk and for mas- 
sive gas accretion streams from the disk to the star enshroud- 
ing the hot coronal plasma. Because the line of sight toward 
DG Tau A is inclin ed against the axis a t a rather small angle of 
37.7 + 3 deg (Eisloffel & Mundlfl998l) . the absorbing mass ac- 
cretion streams must reach high stellar latitudes rather than fall 
toward the equatorial regions. This fully su pports th e commo n 
picture of magnetic accretion (e.g., ICalvet&Gullbrin g i l998). 

Detection of the soft jet component requires low attenua- 
tion by gas, i.e., it is inaccessible in embedded protostars al- 
though such X-rays may locally exist in all jet-driving sources. 
We have indeed found only hard components in our strongly 
obscured targets, and they coincide very closely with the posi- 
tion of the stars. No displaced soft components were identified 
although we re call that a soft jet s ource was reported previously 
for DG Tau A JGudel et all2005t see also Fig.E). 
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